Buckling of semiflexible filaments under compression 



OS 

o 
o 

(N 

X> 

<D 

PL. 
On 



O 



I 

-a 
c 

o 
o 



(N 
> 
On 

o 

(N 

(N 
O 

on 
o 



X 



J. R. Blundell and E.M. Terentjev 

Cavendish Laboratory, University of Cambridge 
JJ Thomson Avenue, Cambridge, CB3 OHE, U.K. 

February 19, 2009 



A model for filament buckling at finite temper- 
atures is presented. Starting from the classi- 
cal worm-like chain model under constant com- 
pression, we use a mean-field approach for fil- 
ament inextensibility to find the complete par- 
tition function. We find that there is a sim- 
ple interpolation formula that describes the free 
energy of chains or filaments as a function of 
end-to-end separation, which spans the whole 
range of filament stiffnesses. Using this formula 
we study the buckling transition of semiflexible 
filaments and find that kinetics plays an impor- 
tant role. We propose that the filament buckling 
is essentially the first order transition governed 
by the kinetics of escaping a local free energy 
minimum. A simple model for the kinetics is 
put forward, which shows the critical buckling 
force for a filament is reduced by a fraction that 
has a universal scaling with temperature with 
an exponent v = 0.56. 

1 Introduction 

There are a number of reasons for the current 
interest in semiflexible filaments under compres- 
sion. Firstly the three structural components 
of the cytoskeleton: microtubules, intermediate 
and actin filaments are all classed as semiflex- 
ible filaments. Experiments on in-vitro micro- 
tubules [U [2] and recent experiments on in-vitro 
dendritic actin [3J, as well as experiments on in- 
vivo actin filaments [I], have shown that fila- 
ments buckling under compression often have a 
crucial role in determining the behavior of the 
cellular network. In addition there is consid- 
erable interest in the mechanical properties of 



carbon nanotubes (CNT's) under compression 
E| , with views to application in drug delivery 
[7J. In light of these considerations, understand- 
ing the behavior of such filaments when subject 
to compressive forces is of considerable impor- 
tance. 

The buckling of rods under longitudinal com- 
pression is a classic problem dating back to Eu- 
ler [S] . For a macroscopic elastic rod with bend- 
ing modulus A and length L that is pinned at 
both ends, Euler derived that above a critical 
force f c = Aw 2 /L 2 the rod can no longer sup- 
port the compressive force and catastrophically 
buckles [HI ED]- The Euler buckling threshold 
only depends on the bending modulus of the 
rod, the compressive (Young) modulus does not 
enter directly into the expression. In fact, a 
simple analysis suggests that the filament ex- 
tensibility can be safely ignored, for homoge- 
neous filaments or hollow tubes alike. Consider 
an elastic tube of outer radius r and the wall 
thickness b (the homogeneous solid rod would 
have b = r) made of a material with the Young 
modulus Y . The bending modulus of such a 
filament is A ~ jYr 3 dr ~ (r 3 b/L)Y [§] and 
the characteristic energy scale of small bending 
is A/L, see Figured! The equivalent strain can 
be achieved by longitudinal compression, with a 
characteristic energy scale ~ (rbL)Y. The ratio 
of bending to compression energies is, therefore, 
(r/L) 2 . For most semiflexible filaments, radii 
are on the order of nm while lengths are on the 
order of microns, so that (r/L) 2 ~ 10" 6 , there- 
fore, most of the imposed deformation will be 
accommodated in bending, so the compressibil- 
ity of filaments can be neglected. We carry on 



1 



here with the simpler problem of incompressible 
filaments. 

If we reduce the scale of rods such that we are 
now dealing with microscopic filaments, there 
comes a regime where the characteristic bend- 
ing energy of the filament A/L is of the same 
order as the thermal energy supplied by the en- 
vironment, ksT. Under such conditions we can 
no longer only consider the mechanical energy 
of the system, but must also include the effect of 
entropy. This is the main purpose of the current 
work. 

The effect of temperature on the buckling tran- 
sition has been considered before, notably in the 
works PH H21 Q31 133!. Odijk evaluates the par- 
tition sum as a semi-classical series expansion 
in three dimensions while Baczynski et. al. [13] 
and Emanuel et. al. [H] consider contributions 
to the partition sum from anharmonic terms in 
two and three dimensions respectively. Hansen 
et. al. [12] evaluate the partition sum by means 
of a saddle point approximation controlled by 
the small parameter l/d where d is the spatial 
dimension and focus largely on the effect of non- 
local interactions along the chain. 
In the present work we take a different ap- 
proach. Invoking a mean-field approximation to 
filament inextensibility, we calculate the com- 
plete partition function and find a simple in- 
terpolation formula for the free energy of sin- 
gle filaments under compression that delineates 
clearly the role of filament bending rigidity and 
the role of entropy, or local thermal motion. 
For arbitrarily small compressive forces we find 
that the free energy function possesses two min- 
ima, one corresponding to the original extended 
filament, somewhat compressed by the applied 
force, while the other describing the completely 
buckled filament. We then propose a simple 
model that describes the first order transition 
from unbuckled to buckled states of the fluc- 
tuating filament. This simple model predicts 
that the critical buckling force is reduced at 
non-zero temperatures from the original Euler 
expression. Our results bear resemblances to 
earlier work [H] in predicting a reduction the 
the critical buckling force at non-zero tempera- 
tures, but our mechanism is different: kinetics 
of overcoming a free energy barrier plays a cru- 



/ +00* 




< (] - x)L >< 


R=xL 



Figure 1: Sketch of a semiflexible filament under 
compression force /, illustrating the notations 
used in the model. 

cial role. 

The paper is organized in the following way. In 
subsection 12.11 we summarize the well explored 
worm-like chain model, which we use with a 
global inextensibility approximation to calcu- 
late the full partition function. We find a simple 
interpolation formula that captures all relevant 
features of this partition function, and the corre- 
sponding free energy of single chains under com- 
pression in a closed analytical expression valid 
across the whole range, from flexible Gaussian 
chains to completely stiff elastic rods; this is 
shown in subsection 12.21 Using this formula we 
examine the effect of compression on the sepa- 
ration of the ends of the chain in subsection 12.31 
and develop a simple model for the kinetics of 
the transition in subsection 12.41 We illustrate 
and discuss the main results of this model in 
section [3j 

2 Theory 

2.1 Semiflexible filaments 

The semiflexible nature of filaments is a result 
of there being an energy penalty associated with 
filament curvature. This energy penalty is con- 
trolled by the bending modulus A with units 
of energy x length. The characteristic bending 
energy of a filament bending over length scales 
of order L is therefore A/L. For a large class 
of microscopic filaments, this energy scale is of 
order &bT, and therefore thermal fluctuations 
cause the filament to bend over length scales 
comparable with L. Such filaments are termed 
semiflexible. In this regime thermal fluctuations 
are large enough to weakly bend the filament 
and thus the filament occupies a middle ground 
between the flexible Gaussian chains (with no 
or weak bending energy) and rigid elastic rods 
whose bending energy far exceeds k^T . 
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The statistical physics of such semiflexible fila- 
ments can be captured by the so called worm- 
like chain model first proposed by Kratky and 
Porod [15] and subsequently developed in fur- 
ther works dH ESI ESI HSl US EEQ • In this model, 
the filament is treated as a space-curve r(s), 
where s parameterizes the contour length along 
the curve. The Hamiltonian is obtained by as- 
sociating an energy that scales like the squared 
local curvature of the space-curve and is propor- 
tional to the bending modulus of the filament A 
which measures resistance to bending: 



H[r(s)} 
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Here the second derivative with respect to the 
arc length s represents the curvature, and we 
also have the constraint that the tangent vector 
has unit length, (dr/ds) 2 = 1 for all s, which 
ensures that the filament is locally inextensible. 
In most situations we are interested in the con- 
strained partition function, which is equivalent 
to the probability that the ends of the filament 
are separated by a span vector R= r(L) — r(0) 
given a filament length L. We use x to denote 
the dimensionless measure of end-to-end separa- 
tion, x = R/L. In the remainder of the paper 
we will be interested in the quantity P(x;L), 
which is the probability of a filament of length 
L adopting a configuration that has its ends sep- 
arated by x. This can be obtained by a func- 
tional integral over the field r(s): 



P(x;L) oc / Vr exp(-H[r]/k B T) 



(2) 



subject to the explicit constraints: 



r(L) -r(0) 



x and 



dj) =1 -< 3 > 



Physically we are summing over the probabili- 
ties of all configurations of the space-curve that 
start at and end at x, and which are locally in- 
extensible. The full expression (j2J) corresponds 
to a non-linear sigma model, which can only be 
solved as a series expansion in the small param- 
eter ksTL/A <C 1, and is therefore only valid 
for very stiff filaments [2CT1 [22]. 



2.2 Mean field inextensibility 

In this paper we relax the constraint of rigid 
local inextensibility to one of global inextensi- 
bility that only counts contributions to the func- 
tional integral whose average tangent vectors 
are of unit length; ( | dr/ds |) = 1. This method 
amounts to a commonly used mean-field ap- 
proach to the problem, as outlined in [T8l [2T] . 
This approach is useful because it makes the 
problem mathematically tractable and produces 
closed forms for the probability distribution, 
while retaining the relevant physical features of 
the filament. This procedure is analogous to the 
transformation from microcanonical to canoni- 
cal ensemble in statistical mechanics. Using this 
approximation the expression for the probabil- 
ity distribution reduces to 



P(x;L) oc J Vr 



-H[r]/k B T 



5[((dr/ds) 2 ) - 1] , (4) 



subject to the boundary conditions (r(L) — 
r(0))/L = x. This functional integral can be 
performed with no further approximation and 
the result can be written as an integral over an 
auxiliary field 4> [21] • Moreover, the probability 
distribution turns out to depend on its variables 
only through a single non-dimensional combina- 
tion of parameters: 
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For a chain in d dimensions the general result 



is: 



P(i) 



sin 



d/2 



(5) 



This expression (for the non-normalized proba- 
bility) uses only the approximation of global in- 
extensibility. It is valid for all values of bending 
stiffness and all extensions of the filament. In 
two dimensions the integral can be solved ana- 
lytically via contour integration [21], and closely 
agrees with previous expressions obtained in the 
literature [20]. In three dimensions the integral 
cannot be performed analytically, however, is 
easily evaluated numerically due to its depen- 
dence only on a single parameter 7. In the 
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Figure 2: Comparison of numerical values for 
the 3-d integral in (J3J (empty circles) with the 
interpolation formula (J5J), plotted as a solid 
line. The interpolation fits very well, which 
is emphasized by the logarithmic axis, in par- 
ticular capturing the correct scaling at large 
and small 7 with the maximum at the correct 
value. For comparison, dashed line shows the 
Ha-Thirumalai expression [18] and the dotted 
line the Marko-Siggia expression |17j . 

past there have been several approximate ana- 
lytical expressions used to represent the proba- 
bility P(x; L) in 3-dimensions, some used quite 
widely [T71 US] . After a new analysis we find 
that the simple interpolation of the form 




P{"i) = exp[— 7r 2 7 



7T7 



(6) 



is a very good approximation to the numerical 
integral (0) for d = 3. This form captures all es- 
sential features of the curve, scaling in the cor- 
rect way for large and small values of 7. A com- 
parison of the interpolation formula ([SD with the 
numerically evaluated values of the integral ([3} 
in 3-d is shown in Figure [21 For comparison, this 
figure also plots the famous and widely used for- 
mulas due to Marko and Siggia [TF] , and Ha and 
Thirumalai [HJ. All these, and other analyti- 
cal expressions agree in the limit 7 — > 0, which 
corresponds to the long or flexible chains, even- 
tually reaching the Gaussian limit. There are, 
however, spectacular disagreements in the limit 
of stiff, or short chains (at 7 > 1), where the 
limit of pure bending elasticity is not present. 
Substituting in the expression for 7, we find that 
within the global inextensibility approximation 
there is a simple expression for the free energy 
F = — kgTlnP(x) of a single chain as a func- 



Figure 3: Non- normalized probability distri- 
bution plots P(x) for filaments of increasing 
stiffness (in order of increasing line width): 
A/k B TL = 0.125, 0.25, 0.5, 1.0, 2.0, 4.0. 



tion of extension, given by 
Avr 2 



F(x) 



2L 



j )+ Mi. (7) 



ttA(1 -x 2 ) 



This form of the chain free energy differs from 
previous results obtained in the literature [TTJ 
IT8"l [TO] in that it clearly delineates the role of 
bending energy and the non-zero temperature 
effects, apart from matching the numerical re- 
sult almost exactly (in fact, we are somewhat 
mystified by the perfect matching of the fit given 
by the expression which suggest that there 
is an underlying analytical route to this integra- 
tion). As one would expect in the retrospect, 
the free energy ([7j) has a term independent of 
temperature that is simply the internal energy 
of the bent filament, and also a term that cap- 
tures the effect of thermal fluctuations. We note 
that this is not simply an entropic term since it 
is proportional to T 2 , it contains a mixture of 
energetic effects and entropic effects at non-zero 
temperature. 

If one takes the limit of a very short filament, or 
equivalently a filament with a very large bend- 
ing modulus A, then the energetic term domi- 
nates over the entropic term, and the free en- 
ergy of the filament is simply the bending en- 
ergy of an elastic rod. In this limit the for- 
mula ([7]) reassuringly recovers the bending en- 
ergy of a macroscopic rod that is derived in the 
small bending regime using a variational princi- 
ple: U « (Atx 2 /L){1 -x) [23]. 
In the opposite limit of small bending modulus 
or a very long filament, the entropic term dom- 
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inates the free energy. In this regime the mini- 
mum of the free energy is at zero extension and, 
expanding about this minimum, in the leading 
order one recovers the free energy of a Gaus- 
sian chain with step length (or the persistence 
length) l p ~ A/ksT. Unlike a true Gaussian 
chain however, higher order terms in our ex- 
pansion would implement the inextensibility of 
the chain that we have imposed. We illustrate 
these limits and the crossover between them by 
plotting the non-normalized probability distri- 
butions for filaments of increasing stiffness in 
Figure [31 

2.3 Filament under compression 

From the compact analytical expression for the 
free energy of a single filament ([6]) we are able to 
write down the free energy for a situation when 
this filament is subject to a constant compres- 
sional force /. In this case there is an addi- 
tional energetic contribution due to the work 
done by the compressive force in moving the 
ends of the filament parallel to the applied force. 
In order that the filament does not rotate, this 
force must be directed parallel to the span vec- 
tor of the filament, and therefore the work done 
is AW —fLx, and so the free energy of a single 
filament takes the form 

rs ,r A ?t 2 , 9, 2(k B T) 2 L 

F(x, f) = fLx + — - 1 - x 2 ) + \ B ' 

1L ttA(1 — x z ) 

(8) 

where the positive / refers to a compressional 
force. Figure H] shows the free energy of a semi- 
flexible filament (at a fixed Ajk^TL = 3) plot- 
ted function of the relative separation for 
increasing compressive forces. 
Examining Figure H] we see there are four 
regimes. For no compression / = 0, the free 
energy possesses a single minimum near full ex- 
tension corresponding to an unbuckled filament. 
This minimum, corresponding to the equilib- 
rium end-to-end distance of a semiflexible fil- 
ament, is at x eq = (1 - 2A; B TL/7r 3 / 2 A) 1 / 2 . The 
linear modulus k (an effective spring constant, 
that is the curvature of the free energy at x eq ) of 
such a filament at small extension / compression 
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Figure 4: F(x; f) plotted for a typical semi- 
flexible filament (A/ksTL = 3) subject to in- 
creasing compressive forces measured in units 
of the Euler critical force f c = Air 2 /L 2 . The 
four regimes of one minimum (unbuckled), two 
minima (unbuckled lower), two minima (buck- 
led lower) and one minimum (buckled) are illus- 
trated by the four curves (low to high) . We stop 
the plots at x = whereas in reality they could 
extend to x < which corresponds to the ends 
"swapping over" and the filament being subject 
to tension. Such effects are not addressed in this 
paper. 

forces about its equilibrium is 

where l p = A/ksT is what one traditionally 
calls the persistence length in longer or more 
flexible chains. Clearly, below a certain bending 
rigidity (or for long enough filaments, L/l p > 
7r 3 / 2 /2) the equilibrium separation x eq = and 
the filament responds to stretching with purely 
entropic elasticity, as classical polymer chains, 
see Figure |3j All of these conclusions are very 
easy to obtain from the new compact analytical 
expression for the free energy (jSJ). 
Returning to semiflexible filaments, for interme- 
diate compressions < / < }\{T) the free en- 
ergy has two minima - one corresponding to an 
unbuckled filament (close to the full extension 
x = 1) and the other corresponding to a buckled 
filament (at x — 0). Provided the compression 
does not exceed f\ , a critical force that depends 
on temperature, the minimum of the extended 
filament is lower in free energy than the buckled 
free energy and so nothing dramatic will hap- 
pen to it: the filament will resist the compres- 
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sion force as an effective spring with a constant 
given by © at / -»■ 0. 

For / > f%(T) there is only one free energy min- 
imum at x = corresponding to the buckled 
state and so the filament cannot sustain the ap- 
plied force. In the region f\{T) < f < /2(F) 
the free energy has two minima, with the mini- 
mum corresponding to a buckled filament lower 
in free energy; the buckling transition may oc- 
cur anywhere in this region depending on the 
relative height of the energy barrier separating 
the metastable extended state. This scenario 
clearly resembles the first-order phase transi- 
tion with a discontinuous jump of the "order 
parameter" and a region of hysteresis. A free 
energy function that possesses two minima of 
equal depth, at / = /i, would be characteris- 
tic of the equilibrium point of the first order 
phase transition. However, as in all discontinu- 
ous transitions, the passage from an unbuckled 
configuration to a buckled one means climbing 
over a free energy barrier. If the system is un- 
able to pass over this barrier in timescales of ob- 
servation, the transition will not occur. The ki- 
netics of this transition will therefore be crucial 
to what is observed experimentally. The force 
fi is not necessarily the force at which buckling 
will occur, instead we are interested in the force, 
which will be called /* where the kinetics of the 
buckling transition changes from being slow to 
fast (in comparison to observation times). 

2.4 Transition Kinetics 

The transition from the extended filament to 
the buckled one means climbing over a free en- 
ergy barrier, the height of which depends on the 
applied compression. The rate of escape over 
this barrier is the classical Kramers escape prob- 
lem [21]. The rate of escape r from a metastable 
local minimum to the true equilibrium state is 
dominated by an exponential term in the barrier 
height AF, i.e. r = r exp(-AF/k B T). If the 
barrier height is far greater than k B T then the 
kinetics of the transition will be very slow. Con- 
versely the escape will take place at an appre- 
ciable rate if the thermal energy is of a similar 
magnitude to the barrier, AF ~ k B T. Strictly, 
this is a delicate kinetic problem involving a hi- 
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Figure 5: A plot of the free energy of a single 
semiflexible filament (A/k B TL = 3) subject to 
a compressive forces about / = /*, focusing on 
the region near the extended state. At / = /* 
(thick line) the barrier is exactly k B T high. For 
compressive forces / < /* the barrier height is 
greater than k B T and the filament in confined 
to the local minimum in separation between x + 
and it remains unbuckled. For forces f > f* 
the barrier height is lower than k B T and we 
assume the thermal energy is enough for the 
filament to "escape" the metastable state and 
buckle. 

erarchy of attempt rates To for different chain 
segments. However, here we take a simple and 
qualitative approach, relying on the strong ef- 
fect of the activation exponential term. Let us 
denote /* the value of external force at which 
AF(/*,T) = k B T, see Figure El At this point 
we would expect there to be enough thermal en- 
ergy in the environment to excite the filament 
out of the local minimum. For compression 
forces above /* we would expect the filament to 
lie in the lower global minimum of the buckled 
state. We therefore propose the rate of escape 
follows an "all or nothing" transition at / = /* 
defined by the conditions 



AF(f,T)>k B T 
AF{f,T)<k B T 



no transition 
transition 



(10) 



where there is equality when / = /*. 
It should be stressed that this "all-or-nothing" 
model is really an oversimplification. The tran- 
sition to a buckled state can proceed even if 
k B T < AF and, equally, a filament can remain 
unbuckled even if k B T > AF, depending on the 
time of observation (or the rate of force appli- 
cation). However the model does capture the 
important physics that it is only in the vicinity 



6 



of AF = kgT that the rate of the transition 
changes from very slow to very fast. More pre- 
cisely, if the rate of escape is of the Kramers' 
form t ~ t exp(— AF/ksT) then the change 
from "slow" to "fast" kinetics occurs at the in- 
flection point of t(T), at AF = 2/c#T; the char- 
acteristic width of this transition is ~ /csT. We 
are justified in assuming the transition rate is 
a step function at AF = ksT if the width of 
the maximum (~ ksT) is small in comparison 
to the overall energy scale of the filament A/L. 
Adopting this simple view of the kinetics there- 
fore places a limit on the range of temperatures 
for which this model will be practically valid. 
The energy scale that determines the height of 
the barrier is A/L. If we are in a regime where 
A/L < ksT then no matter what compressive 
force we impose, the filament is able to explore 
its entire range of extensions in short times. In 
other words, the flexible chain is able to ex- 
plore most of its conformations in the time of 
observation, which is the basis of classical poly- 
mer physics. The filament therefore equilibrates 
rapidly and the assumption that it is trapped 
in a metastable local minimum is a poor one. If 
the system is even more flexible, A/L < 0.5A;£T, 
then it cannot support any compressional force 
at all. Our assumption that the transition is 
a non-equilibrium "all-or-nothing" type transi- 
tion limited by the kinetics of the transition, is 
only valid if A/L > k B T, that is, with no force 
applied, the filament must be confined to the 
unbuckled minimum. For the remainder of the 
paper therefore, we will focus on regime of rel- 
atively stiff filaments where A/ksTL > 1. 

3 Results and discussion 

In this section we discuss the main results of 
our analysis of the filament buckling transition. 
With this "all or nothing" model defined, the 
transition to a buckled state occurs at a critical 
force f*(T) defined by AF(f*,T) = k B T, i.e. 
the force at which the barrier height is ksT. 
This is illustrated in Figure [51 Below /*, the 
filament is confined to the local minimum cen- 
tered on xq of the unbuckled state, and can ex- 
plore the vicinity of the minimum up to points 
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Figure 6: Plots of x+, X- and x® for a semi- 
flexible filament (A/ksTL = 3) subject to in- 
creasing compression measured in units of the 
critical Euler buckling force f c . For increasing 
compression, at / = /* there is no longer a so- 
lution for X- and the filament buckles. On de- 
creasing compression the reverse process occurs 
at a lower force / = /**. Expressions for /* and 
/** are given in the text. 

x + and X- defined as the separations at which 
the free energy is fcgT higher than at the local 
minimum. As the force / — > /*, from below, 
both x + and x_ decrease (see Figure EJ, un- 
til the point / = /* when there is no longer a 
solution for x_ and the filament makes the tran- 
sition to the buckled state. Overall, the range 
between x + and a;_ represents the characteristic 
range of fluctuations of the filament span, which 
is shown as the shaded area in Figure [6j 

Figure [6] shows a plot of and xq for the 

applied force increasing up until the buckling 
transition at / = /*. Above this force the fil- 
ament catastrophically buckles to a separation 
x = and can once again explore the free en- 
ergy minimum up to points that are ksT higher 
in energy than the minimum. If we reverse this 
and now start to reduce the originally applied 
compression force, there comes a point / = /' 
where the minimum at x ~ 1 (unbuckled) be- 
comes lower in free energy than the buckled one 
at x = 0. However, just as for the case of 
buckling, there is now a barrier for departing 
from the metastable buckled configuration (of 
course, assuming that the buckled filament re- 
tains its integrity, which may not be the case 
in great many practical situations). Using the 
same principle as before this transition will oc- 
cur at an appreciable rate only when the barrier 
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Figure 7: A plot of 1 - f*/f c , and inset f**/f c , 
against k B TL/A on log-log axes. The straight 
lines illustrate the simple scaling relations high- 
lighted in the text do indeed hold. 

is of the order of k B T in height. This occurs at 
a different force, which we call /**, that is to 
say the system shows hysteresis (as is typical 
for first-order phase transitions). 
The central result of this work is that we find 
that both f*(T) and f**(T) follow the simple 
scaling laws of the form 

f*/f c = 1 — c {k B TL/A) v 

(11) 

r/fc = d(k B TL/A) r i 

These relations are illustrated in Figure [7] where 
we plot 1 — f*/f c , and inset f**/f c , against 
k B TL/A (which is the only non-dimensional 
combination of model parameters) on log-log 
axes and obtain straight lines. The best fit val- 
ues for the parameters are c = 1.11, v = 0.56, 
d = 0.69, rj = 0.50. As expected, in the limit of 
a very stiff filament (large A) f* — > f c recover- 
ing the macroscopic behavior of a purely elastic 
rod under compression. How the critical buck- 
ling force /* is reduced at finite temperatures is 
the unusual result. We have to emphasize that 
the exponent v = 0.56 is truly universal and re- 
mains independent of what cutoff we choose for 
our "all of nothing" kinetics. Indeed, redefining 
the cutoff to be, say, 2k B T would simply rescale 
our units of temperature, changing the value of 
the constant prefactor c alone. 
For filaments with A/k B TL £s 1 the effects of 
non-zero thermal fluctuations become increas- 
ingly important. The effects described in this 
work should be experimentally observable in mi- 
croscopic filaments such as microtubules, actin 



filaments and carbon nanotubes that are in 
this regime of bending rigidity. For such fil- 
aments this work suggests that the buckling 
force will no longer scale like L~ 2 as it does 
for the Euler buckling of macroscopic rods, but 
like L~ 2 (l — c(L// p ) ' 56 ) where c is a constant 
of order unity and l p = A/k B T is the nominal 
persistence length of the filament. 

4 Conclusions 

We began this paper aiming to examine buck- 
ling in filaments where thermal energy k B T be- 
comes comparable to the energy scales of bend- 
ing A/L. Within the mean field approxima- 
tion of global inextensibility of such semiflexi- 
ble chains we have found that there is a simple 
algebraic expression for the free energy of a fil- 
ament subject to a constant force. This expres- 
sion captures the correct physics of filaments in 
both the flexible limit (flexible Gaussian chains) 
and stiff limit (rigid elastic rods). 
For any non-zero compression, the expression 
obtained for the free energy of single semi- 
flexible chains develops a local minimum cor- 
responding to a buckled state in addition to 
the minimum corresponding to the unbuckled 
state. The buckling transition is therefore de- 
termined by an escape from a local minimum 
in free energy to the global minimum. We have 
adopted a simple "all-or-nothing" kinetic model 
for this transition, in which the transition oc- 
curs if the free energy barrier between the states 
is less than k B T. This suggests that what de- 
termines whether a filament buckles is kinetics: 
how quickly the filament can equilibrate into the 
global free energy minimum. 
Based on this hypothesis, we find that the criti- 
cal buckling force for filaments in a thermal en- 
vironment is no longer the classical macroscopic 
expression f c = ir 2 A/L 2 obtained by Euler. In- 
stead, the critical buckling force is lowered by a 
factor ~ (k B TL/A) ' 56 with a universal scaling 
exponent. This effect should be experimentally 
observable in semiflexible filaments and could 
have an importance for a wide class of filaments 
that are the subject of current research. 
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